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We present a detailed investigation of the structure and activity of extended defects namely 
monoatomic steps on (Ix l)-TiO2(110). Specifically, the two most stable (001) and (111) step 
edges are considered. Employing an automated genetic algorithm that samples a large number of 
candidates for each step edge, more stable, reconstructed structures were found for the (111) step 
edge, while the bulk truncated structures were recovered for the (001) step edge. We demonstrate 
how oxygen vacancies along these defects have lower formation energies than on flat terraces and 
how water and methanol molecules adsorb dissociatively on reduced (111) step edges. Our findings 
are in agreement with earlier experimental results and indicate an important contribution from step 
edges to the reactivity of the TiO2(110) surface. 



I. INTRODUCTION 

Titanium dioxide has attracted much attention due to 
its catalytic and photocatalytic properties, which allow 
a wide range of applications^^. Properties which can be 
extremely different depending on the morphology that 
may range from nanoparticles of various shape and size 
to porous films to polycrystalline materials. All these 
systems have a large fraction of atoms at steps, edges 
and corners, which in the case of metal and semicon- 
ductors are believed to be active sites^^. Although such 
extended defects play an important role also in the physi- 
cal and chemical properties of metal oxide surfaces, much 
less attention has been devoted to these systems^. For 
example, step edges on metal oxides can act as nucleation 
sites for deposited metal clusters^ and are important in 
growth mechanisms and surface reconstructions^JiS. 

Moreover, it is of fundamental importance to under- 
stand the detailed mechanism of the interaction be- 
tween Ti02 surfaces with molecules that are present 
under real conditions (humid air, liquid environment, 
etc.)ii. Water is one of the most studied molecules inter- 
acting with Ti02 surfaces owning to its photocatalytic 
application^ii^. Even in surface science studies that 
use UHV pressures water molecules are always present. 
On the other hand, alcohols interacting with Ti02 sur- 
faces are model systems for the catalytic and photocat- 
alytic oxidation of organic contaminants^^ and methanol, 
being the simplest one, has often been used as probe 
moleculei^ii^. 

Here we study the rutile TiO2(110) surface which is 
considered an important model system for metal oxide 
surfaces -ii^ii^. In particular the unreconstructed (1x1) 
surface is the thermodynamically most stable face of ru- 
tile and is expected to be present on every sample, from 
nanoparticles to polycrystalline aggregates. The activity 
of this surface has been extensively studied both exper- 
imentally and theoretically with focus on surface point 
defects such as bridging oxygen vacancies {Vq)^ hydroxyl 
groups (OH^^) or Ti interstitials in the near surface re- 
gion. However, despite their abundance at the rutile 
(110) surface^ J^, systematic theoretical studies of the 



structure and the activity of extended defects such as 
step edges have never been reported before. 

In this work, we report that reconstructed and reduced 
(ill) step edges on rutile TiO2(110) are mild trapping 
sites for the dissociation of small molecules such as water 
and methanol. Starting from the bulk truncated struc- 
tures of the two most stable (111) and (001) step edges, 
the structures of which are shown in Fig. (TJ we apply 
a systematic search for any other possible structure em- 
ploying an automated genetic algorithm (GA) sampling 
up to 200 different candidates for each step edge. Indeed, 
for the (111) step we find a new configuration that is more 
stable with respect to the bulk truncated one. When the 
new structure is considered, we find that oxygen vacan- 
cies along step edges have a low formation energy, mean- 
ing that this step is easily reducible. We further demon- 
strate how water and methanol dissociate on reduced and 
reconstructed (111) steps with binding energies of about 
0.5 eV lower with respect to the highly reactive bridging 
Vo' Our findings are in agreement with published wa- 
ter temperature programmed desorption (TPD) experi- 
ments and we propose a new assignment for methanol 
TPD measurements. 



II. METHODS 
A. Computational details 

The calculations were performed at the Generalized 
Gradient Approximation (GGA) level of DFT using 
the Perdew-Burke-Ernzerhof (PBE) exchange correlation 
functional^^ as implemented in the Grid-based Projector 
Augmented- Wave (GPAW) program^^'^^. Electrons are 
described using the projector augmented wave (PAW) 
method in the frozen core approximation^^. 

The (001) and (111) steps have been modeled with vici- 
nal rutile (230) and (451) surfaces, respectively. Periodic 
boundary conditions have been applied only along the 
surface plane (x,?/), whereas a vacuum of at least 7 A 
between the surfaces and the walls of the cell was set. 
A grid spacing of no more than 0.20 A was used. The 
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area of the vicinal surface times the number of k-points in 
the Brihouin zone is always bigger than 640 A^. A four 
electron setup (4e) has been used for Ti atoms, which 
consider the 3d and 4^ electrons in the valence shell. The 
equilibrium lattice parameters of the bulk rutile phase 
are a = 4.71 A and c = 3.01 A. 

The graph in Fig. [2] has been calculated on a 
Ti37074(451) slab. We have checked the accuracy of our 
calculations with the 4e~ setup by employing a 12 elec- 
tron setup (12e~), which includes also the 35 and the 3p 
electrons in the valence shell. The equilibrium lattice pa- 
rameters calculated for the 12e~ setup are a = 4.65 A and 
c = 2.96 A. The results using the 4e~ and 12e~ setups 
shows good correspondence (c/. caption of Fig. [2]). 

Oxygen vacancies have been created starting from a 
stoichiometric Tiii8 0236(451) slab, which corresponds to 
a (2x1) unit cell. A (451)-(lxl) unit cell has been used 
to simulate a high concentration of oxygen vacancies. 
Such configuration results in a step Vo every O row or 
one Vo per ~7 A. 

Adsorption of water and methanol on a TiO2(110) ter- 
race has been preformed on a four Ti02 trilayer slab using 
a (4x2) unit cell. The slab is asymmetric with the adsor- 
bate on one side only. The formation energy of a bridging 
Vo with respect to an O2 molecule in the gas phase is 
3.26 eV. The adsorption of water and methanol on a step 
oxygen vacancy has been modeled in a Ti590ii7(451) 
slab. 



B. Genetic Algorithm 

When investigating a possible step reconstruction the 
atoms in the two rows of Ti02 closest to the step are 
considered the target atoms for an optimization with a 
newly implemented genetic algorithm (GA)^^. The GA 
starts from a population of 15 randomly generated candi- 
dates. A fixed population size of the 15 best candidates is 
maintained throughout. To preserve the stoichiometry in 
the pairing processes each atomic species is paired sepa- 
rately with the cut-and-splice operator Mutations are 
introduced to ensure a diverse population with the pri- 
mary mutation being a permutation of half of the atomic 
numbers in 30% of the newly generated candidates. A 
similarity criterion based both on energy differences and 
a pair distribution function is used to limit the num- 
ber of duplicates in the population. It is however not 
possible to fully hinder duplicates and the GA search is 
considered converged when small variations of the best 
candidate start dominating the population. In order to 
sample a large number of structures we have used a faster 
approach during the DFT optimization step in the GA. 
This method uses a Linear Combination of Atomic Or- 
bitals (LCAO) to describe the wave functions^^. 




FIG. 1. (a) Top and (b) tilted side view of the TiO2(110) 
surface with bulk truncated (001) and (111) step edges. 



III. RESULTS AND DISCUSSION 
A. Step Structures 

Figure [T] illustrates the rutile TiO2(110) unrecon- 
structed surface with bulk truncated steps running along 
the [001] and [ill] directions. Indeed, typical STM im- 
ages show that terraces are bound by monoatomic (001) 
and (111) steps of about 3.2 A height Steps that run 
along other directions are not observed under normal 
condition^. In the present work, we have performed 
individual searches for the optimum structures of steps 
along the [111] and [001] directions. We first report the 
result for the (001) step for which we find the most stable 
step structure to be the bulk truncated one. 

(001) step edges: Unreconstructed (001) step edges can 
be created with two different terminations: Obr rows and 



Ti troughs, (001) o and (001) n respectively, see Fig. [H 
Based on STM images Diebold et al}^ have proposed 
that (001) steps are O^r terminated. We have apphed to 
the (001) steps our GA search in order to find more stable 
structures than the bulk truncated ones. However, after 
sampling up to 140 configurations, the two bulk trun- 
cated structures are recovered as the most stable ones. 
From our calculations the (001 )o step is more stable by 
less than 0.05 eV per unit cell. However, this energy dif- 
ference is at the limit of our precision and therefore we 
are not able to confirm the experimental observation of 
(001) o step to be the most stable one. 

(Ill) step edges: We now move to (111) step edges that 
are expected to be present in a higher concentration^^. 
In Fig. [1^ the (451) surface unit cell used in the search 
for the most stable (111) step is shown. The cell offers 
the smallest possible periodicity along the step direction, 
and separates the steps by about 13 A via fiat TiO2(110) 
terraces containing four TiO2(110) unit cells in each row. 
A view at the (451) unit cell in Fig. [1^ reveals that the 
step edge contains two inequivalent Ti atoms, five- and 
fourfold (5/- and 4/-) coordinated (Fig. [1)3), meaning that 
in order to add one row to the step and recover the same 
step structure two Ti02 units must be added. This in 
turn means that the most stable step may either result 
from steps having an odd or an even number of such 
added Ti02 units. 

The growth of the (111) step is illustrated in Fig. [2l 
Starting from the bulk truncated step edge. Fig. [2^, one 
unit of Ti02 is adsorbed on the terrace in the structure 
reported by Park et al?^^ and Wendt et al?^ , Fig. [2)3. 
With respect to bulk Ti02 hosting the Ti02 on the ter- 
race is unfavorable by about 1.3 eV (Ei = 1.32 eV). Mov- 
ing the Ti02 unit to be part of the (111) step and after 
the search with the GA involving 194 structural candi- 
dates, the structure of Fig. [2]: resulted as the most stable. 
Here, the Ti02 unit is incorporated into the step edge 
and 1.65 eV is gained rendering this configuration the 
thermodynamically preferred state of the (111) step edge 
under condition of stoichiometry by 0.33 eV (E^). We 
shall refer to this structure as (lll)^ and return to the 
question of non-stoichiometric conditions in Section llHBI 
One characteristic of the (lll)i? structure is that the Ti 
atom added to the step edge (purple in Fig. [St) appears 
at a height, 2.16 A, which is intermediate between that 
of the two terraces forming the edge, = A and 3.26 A. 
In Fig. [3^1 we report the energy of all the structures sam- 
pled by the GA. The six most stable ones are illustrated 
in Fig. [3)3. 

To complete the growth cycle and prove the consis- 
tency of our approach, a 2"^ Ti02 unit was added on the 
terrace. Fig. 01, which costs about 1.4 eV (E2 = 1.39 eV). 
Moving this Ti02 species (cyan color in Fig. [2j) to the 
step edge and employing our search with the GA covering 
more than 100 structural candidates, the bulk truncated 
edge structure. Fig. [2^, was recovered. The energy of the 
structures sampled by our GA within 10 eV are reported 
in Fig. [St, and the most stable structures are shown in 
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FIG. 2. Growth of a (111) step edge when going from (a) 
a bulk truncated (111) step to the same structure but with 
two extra Ti02 units. Configuration (c) represents the most 
stable structure, (11 1)^. (e) Potential energy diagram for 
the six structures. The energy values calculated with a 12e~ 
setup are: Ei=1.49 eV; Eh=-0.19 eV; E2=1.40 eV. (f-g) Top 
view of the (111) and (111)^ step edges. The coordination 
numbers of the Ti and the O atoms are indicated counting 
Ti-0 bonds of length up to 2.5 A. 



Fig. [an. 

We emphasize that the two structures. Fig. and[2t, 
have resulted directly from the GA search without any 
constraints or initial guesses of the configurations {of. 
Fig. [3]), confirming the high stability of such step edges 
and proving the reliability of our approach. 

We now address the reasons for the particular high sta- 
bility of the (lll)i? step with respect to the (111) step. 
The two structures for the steps running along the [111] 
direction are shown enlarged in Fig[2]F-g. From the struc- 
ture reported in Fig. [2^ we note that the two oxygen 
atoms of the extra Ti02 unit, are in their bulk positions 
and are twofold coordinated. The extra Ti atom is six- 
fold coordinated in an octahedral site at half height be- 
tween the upper and the lower terrace. Furthermore, the 
new in-plane O atom on the upper terrace establishes the 
terrace-like fivefold coordination of the 4/-Ti atom along 
the bulk truncated (111) step. The new {111)r step that 
we have discovered thus optimizes the coordination of all 
Ti atoms resulting in a very stable structure. 
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FIG. 3. Energy distribution for all the configurations sampled by the GA for the (a) (lil)i? and (c) (111) step edges. The zero 
in energy corresponds to the bulk truncated (111) structure. The size of the dots is proportional to the number of times the 
configuration has been successfully discovered. (b,d) Ball models of the best six structures for the (111)^ and the (111) step 
edges. The Ti02 units included in the GA are in highlighted colors. Energy values calculated with the LCAO (red) and the 
grid-based (black) methods are reported. 



B. Reduced Step Edges 

Having establish a new structure for the (111) steps we 
now turn to the study of oxygen vacancies (Vq) formed 
along such step edges. A reduced Ti02 surface typi- 
cally contains bridging Vo in a concentration from 2 to 
10% ML which are easily imaged by STM as bright pro- 
trusions along the dark Obr rows. The presence of Vo 
at step edges with different coordination with respect to 
bridging Vo is intriguing and a different activity can be 
expected. Indeed, we find that creating a Vo along this 
step is easier compared to the creation of a Vo on the 
terrace. 

In Fig. [4] the formation energies for a series of Vo cre- 
ated by removing different O atoms in our model of the 
(lll)i^ step are reported. Oxygen atoms numbered 2, 
3 or 4 are representative of terrace Obr atoms. A Vo 
created by removing one of these atoms costs about 3.1- 
3.3 eV. The oxygen atoms that are marked from 5 to 9 
are instead characteristic of the (lll)i^ step. The results 
show how steps are clearly highly reducible with Vo at 
position 5 having the lowest formation energy, 2.03 eV. 
Also removing the O atom number 6 costs the same 
amount of energy but in this case the structure undergoes 
a strong relaxation and a Vo in position 5 is created. We 
explain the low formation energy of a Vo in position 5 
in terms of electronic states associated with the O atoms 
prior to the Vo formation. In Fig. the projected den- 



sity of states of the 2p orbitals of the O atoms of interest 
are shown. It is striking that the two O atoms that are 
easy to remove are the ones that contribute to the top of 
the valence band with a well defined state of 2p character. 
Fig. Hh. This state is localized on atoms number 5 and 6 
as depicted in Fig. HJ). Moreover, it has an anti-bonding 
character with a nodal plane between the two atoms, ex- 
plaining the lower energy cost of removing any of them. 
The quoted Vo formation energies apply to the removal 
of every 2"^ atom of the given type along the step edge. 
When increasing the coverage and removing one O atom 
in every O^r row the Vo formation energies remain low 
(2.25 eV for site no. 5) indicating that (lll)i? steps can 
accommodate a high density of oxygen vacancies, one Vo 
per ~7 A. 

Based on these results, we conclude that a reduced 
TiO2(110) sample that shows even a low concentration 
of bridging Vo at the terraces, must accommodate a high 
number of Vo along {111)r step edges. 



C. Adsorption of Water and Methanol 

Given the structures for reduced steps running along 
the [111] direction we have tested the activity of such 
extended defects towards the adsorption of water and 
methanol and compared to the activity of the well known 
reduced and flat (110) terrace. Fig. [5l Both water and 
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FIG. 4. (a) step edge. The O atoms where Vo are 

created are marked in yellow, cyan and purple, (b) Highest 
occupied Kohn-Sham orbital. The isodensity value displayed 
is ±0.4 A~^/^. (c) Total density of states and 2p projected 
density of states of the 9 O atoms marked in (a). On the right 
side of the plot the corresponding Vo formation energies are 
reported. 



methanol molecules dissociate and bind strongly at bridg- 
ing Vo (Fig. El lower part) in agreement with previous 
theoretical and experimental worksi^*^— . In contrast 
to a molecular adsorption along Ti troughs characterized 
by lower binding energies, bridging Vo are considered the 
only trapping site where these molecules can dissociate. 
In Fig. [5] the most stable configurations for the adsorp- 
tion of H2O and CH3OH into a step Vo are illustrated. 
Both molecules dissociate with binding energies of 0.91 
and 1.02 eV for water and methanol, respectively. 

The results for dissociation of water at step edges are in 
very good agreement with earlier water TPD experiments 
that show a recombinative desorption peak at tempera- 
tures lower than the one at vacancy sites^*^. However, 
the presence of this peak in H2O TPD seems to be de- 
pendent on the experimental condition (e.g. see TPD 
spectra in Ref j^^i^^ compared to the one in Ref.^^ where 
this intermediate peak is not discussed). In addition, 
hydroxyl groups at step edges resulted from water disso- 
ciation have been reported to be observed by Suzuki et 

Our results suggest a similar behavior of water and 
methanol. We predict methanol to dissociate at re- 
duced (111) it: steps. This result is compatible with 
methanol TPD experiments performed by Henderson et 
al.^^ . The TPD spectra for methanol adsorbed on a 
TiO2(110) surface shows three distinct states at 295, 350 
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FIG. 5. Dissociative adsorption for a water and a methanol 
molecule in an oxygen vacancy at step (top) and in a bridging 
oxygen vacancy (bottom). 



and 480 K. The low and the high temperature peaks 
are assigned to molecular methanol adsorbed along Ti 
troughs (E5=0.84 eV) and to dissociative methanol ad- 
sorbed at bridging Vo sites (E5=1.51 eV, Fig. [5]), re- 
spectively. We propose the methanol dissociated at va- 
cancy sites along step edges, the configuration shown in 
Fig. O to be responsible of the third peak at 350 K 
(E5=1.02 eV). This assignment is compatible with the 
appearance of this peak at low coverage and with the au- 
thors' conclusion that methanol desorbing at this tem- 
perature is formed from recombination of dissociated 
methoxy species. 

IV. CONCLUSIONS 

In summary, we have provided a clear understanding 
of the local structure of step edges on rutile TiO2(110) 
surface. An automated genetic algorithm was employed 
leading to the finding of a new structure for the (111) step 
edge that differed from the bulk truncated one. Know- 
ing the correct structure of step edges is of fundamen- 
tal importance and can have many consequences for sev- 
eral reasons. Here, we have shown that oxygen vacancies 
can be extremely stable at step edges. A reduced ru- 
tile sample thus hosts not only point defects, like oxygen 
vacancies or Ti interstitial, but also highly reduced step 
edges. This conclusion can easily be extended to other 
reducible metal oxide system and can strongly influence 
the activity of the near step region. For example, wa- 
ter and methanol molecules dissociate at reduced (111) 
step edges as in bridging Vb, but with lower binding en- 
ergies. This implies that while point defects at terraces 
become polluted, point defects at step may engage in a 
catalytic cycle. The results for water are in good agree- 
ment with earlier experimental results where it has been 
speculated that step edges may play a role in the wa- 
ter dissociation^ii^i^. For methanol adsorption we rein- 
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terpret previous experimental results^^i^ and propose a 
similar behavior for higher alcohols. 

Our findings also suggest the possibility of tuning 
the activity of the rutile TiO2(110) surface allowing to 
decrease the temperature needed for the desorption of 
strongly bounded species. For example, changing the 
experimental conditions or the surrounding environment 
can result in a different morphology of the surface^^'^^ 
and in particular in a different ratio between flat (110) 
terraces and step edges. Moreover, vicinal surfaces that 



expose a high concentration of step edges oriented along 
specific directions can be employed^. 
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